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Abstract. Let / (z, z) be a positive bi-homogeneous hermitian form on C", of de- 
gree to. A theorem proved by Quillen and rediscovered by Catlin and D'Angclo states 
£^ , that for N large enough, (z, z) N f (z, z) can be written as the sum of squares of homo- 

^ ' geneous polynomials of degree m+N. We show this works for N > Cf((n+m) logn) 3 

where C/ has a natural expression in terms of coefficients of /, inversely proportional 
to the minimum of / on the sphere. The proof uses a semiclassical point of view on 
which l/N plays a role of the small parameter h. 
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Let / = f(z, z) be a bi-homogeneous form of degree m > 1 on C": 



f(z,z):= J2 c ^z a z^ zeC n , c aP EC. (l.r 

00 ' \ a \=\p\=m 

Here n > 2, a = (ai,...,a n ) e N n , \a\ := a x + ... + a n , z a := z^...z° n . The following 



theorem was proved by Quillen in 1968 [9], and rediscovered by Catlin and D'Angelo 
in 1996 [2]: 

Theorem 1. Suppose f is given by (1.1) and that 

f(z,z)>0, z^O. 
Then there exists N such that for N > N 

\\z\rf(z,z) = Y.\ P i N ( z )\ 2 > N 2 :=£M 2 ' (1-2) 

j=l j=l 

where P^{z) are homogeneous polynomials of degree m + N, and d^ = (™+™ +A ^ i s tke 
dimension of the space of homogeneous polynomials of degree m + N. 

This result can be considered as the complex variables analogue of Hilbert's 17th 
problem: given a multivariate polynomial that takes only non-negative values over the 
reals, can it be represented as a sum of squares of rational functions? The positive 
answer to this original question was given by Artin in 1926 [1]. For a survey of recent 
work on the hermitian case see the review paper by DAngelo [3]. 

In this paper we give the following quantitative version of Theorem 1: 

l 
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Theorem 2. Let f satisfy the assumptions of Theorem 1 and define 

1/2 

(1.3) 




Then there exists a universal constant C such that (1.2) holds for 

N > C ^-(m + n) 3 log 3 n. (1.4) 

The proofs of Quillen [9] and Catlin-D'Angelo [2] are based on functional analytic 
methods related to the study of Toeplitz operators. The existence of iVo such that 
(1.2) is satisfied is obtained by a non-constructive Fredholm compactness argument - 
see [7, Section 10] for outlines and comparisons of the two proofs, and also [4] for an 
elementary introduction to the subject. 

Here we take a point of view based on the semiclassical study of Toeplitz operators 
- see [11, Chapter 13] and references given there. Our proof of Theorem 2 is a quanti- 
tative version of the proof of Theorem 1 given in [11, Section 13.5.4]: the compactness 
argument is replaced by an asymptotic argument with N — 1/h, where h is the semi- 
classical parameter. The symbol calculus for Toeplitz operators allows estimates in 
terms of h which then translate into a bound on N. 



Better bounds on N obtained using purely algebraic methods already exist and it is 
an interesting question if such bounds can be obtained using semiclassical methods. 

In the diagonal (real) = if a 7^ /?, Theorem 1 is equivalent to a classical 

theorem of Polya - see [7, Section 10.1]. In that case a sharp bound on N was given 
by Powers and Resnick [6]: 

m(m-l)A(f) ~ , fa!, ,1 , , 

2 A (/) \a\=m [ml J 

It is remarkable that the bound does not depend on the dimension n. To compare 
this bound to the bound obtained using semiclassical methods, we note that in the 
diagonal case, the spectral radius used in Lemma (3.1) is given by A(/). Hence an 
easy modification of that lemma leads to the bound 

N> Y^(n + m) 3 log 3 n, (1.6) 

which is weaker than the bound (1.5) from [6], roughly by a factor of m(l + n/m) 3 . 

In the complex case, To-Yeung [10, Theorem 1] gave an algebraic proof of a better 
bound than the one provided by our method in Theorem 2. They show that 

y\ttf f) 

N > nm{2m - 1) \ - n - m, A«(/) := sup \f{z, z)\. 
l0 g^ A UJ |2|=i 
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The common feature of all these bounds is the denominator A(/) and the standard 
example of \zi\ 4 + |^ 2 | 4 ~ c| ^! | 2 1 ^ 2 1 2 , < c < 2, z G C 2 (see for instance [11, Section 
13.5.4]) shows that the 1/A(/) behaviour is optimal. 

In Putinar's generalization of Polya's theorem [8], a much larger bound was given 
by Nie and Schweighofer [5]: 



for some c > 0. 

The paper is organized as follows. In Section 2 we recall various basic facts about 
the Bargmann-Fock space and Toeplitz quantization. Section 3 presents the basic 
inequality which leads to a bound on N. Section 4 provides quantitative estimates 
on the localization of homogeneous polynomials in the Bargmann-Fock space, with a 
stationary phase argument given in the appendix. The proof of Theorem 2 is then 
given in Section 5. 

Notation. We denote (x,y) for G C™ the euclidean quadratic form on C n (not 
the hermitian scalar product): (x,y) := Y^7=i x iVi- ^ OT z = ( Zl -> ■■■■> z n) £ C™ we define 
as the standard hermitian norm: ||z|| 2 := *YTi=\ z i z ~i = ( z > z )- The measure dm(z) 
denotes the 2n- dimensional Lebesgue measure on C n . The space of homogeneous 
polynomials of degree M is denoted Vm- Finally, for two quantities A, B, we write 
A > B, if there exist a (large, universal) constant C, such that A > CB. 

Acknowledgments. We would like to thank Mihai Putinar for encouraging us to 
pursue this project and John D'Angelo for providing useful comments and and refer- 
ences, especially reference [10]. The partial support by the National Science Foundation 
for partial support under the grant DMS-1201417 is also gratefully acknowledged. 

2. Preliminaries: Bargmann-Fock space and Toeplitz quantization 

Quillen's original proof of Theorem 1 used the Bargmann-Fock space - see [7, Section 
10], [9] and [11, Section 13.5.4]. We modify it slightly by introducing a semiclassical 
parameter h and considering the subpace of homogeneous polynomials of degree M, 



We denote the resulting space by L\. The closed subspace of holomorphic functions is 
denoted by H$. The measure exp(— \\z\\ 2 /h)dm(z) will sometimes be written as dG(z). 




(1.7) 



A Hilbert space Bargmann-Fock norm on Vm is given by 




and we can extend this norm to any function u such that 
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The Bergman projector 11$, is the orthogonal projector L% — > H$> and to compute it 
we identify an orthonormal basis of H§. The following standard lemma is a rephrasing 
of [11, Theorem 13.16]: 

Lemma 2.1. Let us define 

i / i y /2 

r^-T (2-1) 



{irh) n/2 \h^a\ 
Then 

(i) The set of f a 's is an othonormal basis on Hq,. 
(ii) The Bergman projector 11$ can be written 



where 



ll$ii (z) = II (z, w) u (w) dm (w) 
R{z,w) : = -^^exp f~ ((z,w) - \w\ 2 ) 



To connect the study of positive bi-homogeneous forms to Bargmann-Fock space, 
we recall the standard result (see [11, Lemma 13.17]): 

Lemma 2.2. A bi-homogeneous form of degree m can be written as a sum of squares 
of homogeneous polynomials, 



k 

f(z,z) = Y,\PAz)\ 2 , PM= E^> 

J=l |a|=m 

if and only if the matrix {^ap)\ a \ = \M- m ^ s positive semidefinite. 

Thus to prove Theorem 1 we need to show that the matrix of the hermitian form 
(z, z) N f (z, z) is positive for N large enough. Let us compute this matrix. Since 



(Z,Z) N _ ^ z"*" 

I . . I AT I 



(z,z) N f(z,z)= C -^-z^z^= <p^> 



where 



U! 

\a\ = \P\=m \j\ = \p\ =m +N 

ImI— n 



E % \ P \ = h\=N + m. (2.2) 

a+p,=p V" 



f-n 



The following essential idea comes from the work of Quillen in [9]. It relates the 
positivity of the matrix (2.2) to the positivity of a differential operator. 
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Let Pf be the following differential operator 

P f= Yl c a pz a (hd z f : ^ H*. (2.3) 

\a\ = \f3\=m 

Since / is real, c^g = cp a . Thus the formula (2.5) shows that Pf is self adjoint. Let us 
explain now how the positivity condition and the operator Pf are related. 

A simple calculation (see [11, Section 13.5.5]) based on the definition and (2.4) shows 
that for all u, v G V m+ N, 

(P f u, v) Vm+N = ^h n+N+2m £ p\l\c%ujv p 

| 7 | = |p|=m+AT 

where u p G C, u 7 G C, are given in 

u = u p z p , v = v 1 z 1 . 

\p\=m+N \-y\=m+N 

Thus proving that the matrix (2.2) is positive definite is equivalent to proving that Pf 
is a positive operator on V m+ N- To make this quantitative we use the following lemma 
which is an application of a more general formula given in [11, Theorem 13.10]: 

Lemma 2.3. Let 11$ be the orthogonal projector from L\ to Hq,. Then 

Pf\p m+N = E co^n. (^-) . (2-4) 

\ct\=\P\=m 

and 

Pf\ Pm+N = Tl m {z,z)^ (2.5) 



where 



^ hi 

= £-H A ) : '/(*,*)■ (2-6) 
i=o J ' 



Using (2.5), positivity of P/ on Vn+tti. follows from inequality 



(Yl^qll^u, u) Vm+N > c\\u\\ 2 L 2 , ueV 



* 



N+mi 



for some constant c > 0. But since n$tt = u and I1J = 11$, it suffices to prove that 
for all u G Vn+iu, with L|-norm equal to 1, 

(q(z,z)u,u) L 2 9 > c, (2.7) 
and (2.7) is the starting point of our work. 
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3. The basic estimate 

We define the ring Q e as 

tt £ := {z e C n , 1 - e < \\z\\ 2 < 1 + e}. 
For u G Vn+tu. with L|-norm equal to 1 we have 

(qu, u) L 2 = q(z, z) \u (z) |V l|z||2//l dm (z) 



Q 



(z,z)\u(z)\ 2 e-^ 2/h dm(z)+ / q (z, z) \u (z) \ 2 e-^ 2/h dm (z) 



> min q \\u\\ L 2 - u L m ^ 



q{z,z) \u(z) | 2 e-W /*dm(z) 



C"\fi e 

Recalling (2.6) we see that 

/■ (-|A) J /( 



j=0 



2;, z 



in 



> 



max 



z ||2(m-i) 
1 



z J n 2 e 



j! n * V Ikll 2 ^ - ^ 



(!a) j 7(, 



\H u \\Ll(n £ ) 



h? 



>-^2-^E e (h,m + N,m-j)m^ (±A) J f ( 



-2, Z 



3=0 

where the quantity E s (h, M, k) is defined as 

E e (h,M,k):= sup 

u&Vm,\W\\v m =1 

and where we used the homogeneity of A J '/ of degree 2 (m — j). 
Rearranging the terms we obtain 



1*11 «IU|(n e )> 



(qu, u) L 2 > I (1 — E e (h, m + iV, 0)) min g 

w 



- Yl (h,m + N,m- j) max (\ A) J / ( 



Moreover, 



mm q (z, z) 

l-2e<||2|| 2 <l+2e 



3=0 



min yiLri&yffz 

l-2e<\\z\\*<l+2ef^j\ ^ 4 > J V 



>(l-2,) m A(/)-^ T (l + 2e) 



max 



(3.1) 



(3.2) 



(3.3) 



(JA)'/(M) 



A QUANTITATIVE VERSION OF C ATLIN-D ' ANGELO-QUILLEN THEOREM 



7 



Z,Z, 



We see that we need an upper bound for max QA) J / ( 
following 

Lemma 3.1. We have the estimate: 

rnax|(iA) j /(^^)| < (nm 2 ) J A(/), 

where A(/) is defined in (1.3). 



and that is given in the 



(3.4) 



To explain the proof we note that since / is a bihomogeneous form of degree m, 
A k f is a bihomogeneous form of degree m — k. If we have estimates on /, and if 
we find an explicit relation between estimates on / and A/, related to the bound on 
max|| 2 || =1 |/ (z, z) |, a recursion procedure will give (3.4) 

Proof. For z G C n satisfying ||z|| = 1, put z = (rie ldl , ...,r n e t9n ), with J^rf = 1. Then 



\f( Z ,z)\< J2 

\a\=\/3\=m 



<m\ „, mi 



\a\=\P\=m 

(CR,R), C 



c a p\\ —r 
ml V ct\ 




ml 



\Caf3 



R :-- 



Since 



E 

|a|=m 



ml 



:2a 



\a\ = \P\=m 



I 2m 




\a\=m 



ft! 



we have 



\f(z,z)\<{CR,R)<p(C), 



where p(C) is the spectral radius of C. The spectral radius can be estimated by A(/) 
given in (1.3): we write C = UDU" 1 , where D and U are diagonal and orthogonal 
matrix, respectively. Then 

A (f) 2 = tr (CC*) = tr (UD 2 ^ 1 ) = tr (D 2 ) > p(C) 2 



and hence 



max\f{z,z) \ <A(f) 



2 =1 



(3.5) 



We now need to find a relation between A (/) and A (|A/J. Let D := (d lp ) be the 
matrix of the bi-homogeneous form |A/, and let us chose 7, p with |7| = \p\ = m — 1. 
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Denoting by c a p the entries of C we obtain 

1 d 1 d p n d d n 

^ = E a?^/ (°> °) = E (7* + 1) (P* + 1) <W*, (3.6) 

i=l i=0 
n I 

= E (7< + 1) (Pi + 1) /7 '., =^fC J+ehP+ei (3.7) 

£f V(7 + ei)!(p + ei)! 

= -pri E v 7 ^ + 1) (Pi + l)5 7+e . p+ej (3.8) 
(m — 1)! 2 _ 



i=i 



If we put d 7P := y/"f\p\d lp / (m — 1)!, and denote the correspoding matrix by £), then 

n 

dyp < m 2 C 7+eil p+e a , (3.10) 



i=l 



and 



x 2 



A (i A /) 2 = E <^ E ™ 4 E^+J (3.H) 

|7| = |p|=m— 1 | 7 | = |p|=m— 1 \ i=l / 

n 

<m 4 n E^,p^, (3-12) 

|7| = |p|=m— 1 i=l 

< nm 4 ■ nA (f) 2 . (3.13) 
An easy recursion then gives 

A((iA)V)<(nm 2 ) J A(/). 
and inequality (3.5) applied to QA) J / instead of / proves the lemma. □ 
The lemma and the lower bound stated after the inequality (3.3) imply 

mm q(z,z) > \{f) {I - 2e) m - K{f)^- {nm 2 h) j {I + 2e) m " ] . (3.14) 
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This combined with (3.3) leads to the basic inequality: 
{qu,u) L 2 > 



;i _ e £ (h, m + N, 0)) f A (/) (1 - 2e) m -A(/)fi (nm 2 h) j (1 + 2e) m ' 

7=1 jl J (3.15) 



m 

A (/) - {nm 2 h) j E £ (h, m + N, m — j) 



3=0 J ' 

All the work that follows is aimed at finding ho such that for h < ho the right hand 
side of (3.15) is positive. 

4. Estimates on E £ 

Our goal in this section is to prove that the quantity E £ (h, M, m) roughly decreases 
like exp (— Me 2 ), under some assumptions relating e,h,M,m,n. It is essentially due 
to the fact that the homogeneous polynomials are localised in L|-norm around the 
sphere S 2n ~ 1 C C n , with l/h ~ M - see [11, Theorem 13.16, (ii)] for an explanation 
of this using the harmonic oscillator. Here we prove 



Lemma 4.1. Let e, h, m, n, M > and let us call 

a := h(M + m + n- 1). (4.1) 

Assume that 

->a>l, I>e>4(<7-1). (4.2) 



Then for E e defined in (3.2) we have 



1 / Me 



E £ (h, M, m) < h m (M + m + n) 2n+m - exp ) • (4-3) 

Proof. Let U^, 1 be the projection from L| to Vm- For u £ Vm- To estimate the right 
hand side in (3.2) we note that 

iiikir«iiU^) = <«' n »ii^ii ami o. II » u >^ 

rAfii „l|2mn ttMii IL.II2 



^ IITTJ" II TTiVJ II ||„, 

< 1)11$ in e ll$ ■ \\u 



Hence it suffices to estimate the norm operator ||n|f ||z|| 2m lfi s n|f ||, and for that we 
will use the following standard variant of Schur's Lemma: 

Lemma 4.2. Let (X, fi) be a measure space, K : L 2 (X) — > L 2 (X) a selfadjoint operator 
with kernel k, that is 

Ku{x) = / k(x,y)u(y)dfi(y). 
J x 

Assume that there exists an almost everywhere positive function p on X and A > 
such that 

\Hx,y)\p(y)dfi(y) < Xp(x). (4.4) 

x 
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Then \\K\\ < X. 

To apply the lemma we first construct the kernel of the projector 11^ = Y^\ a \=M /<*/a 
where f a was defined in (2.1), and /* is the linear form (f a , Writing 

U%u(z):= I Tl M (z,w)u(w)e-^ /h dr 

we have 



\a\=M 



V —- — — z a vT = - V -z a 



1 V- l^a M 



M 



(7rh) n \h M a\J n n h n+M ^ a\ M\n n h n+M ' 

\a\=M v ' v 7 \a\=M 

It follows that the integral kernel of K = H$ WzW 2171 !^^!!^ 1 with respect to the Gaussian 
measure dG(z) := exp(— \\z\\ 2 /h)dm(z), k, is given by 

kM = L, m^m^ U^O- 

This suggests natural choice of the weight p = \\z\\ M in lemma 4.2, and we need to 
estimate the corresponding parameter A in (4.4). For that, we need an upper bound 
on the integral 

\k(z,w)\\\w\\ M dG(w). 

c 

An application of the Cauchy-Schwarz inequality inequality gives 

c ,*(,, „)i W r««j W < 1 1 ikir S^^l^ Kir^odGM 

Thus it is sufficient to estimate the following integrals: 

\2M r ||/M|2M+2m 



A polar coordinates change of variables, followed by a substitution t = r 2 /h, gives 

(4.6) 



10271—11 poo I c2n-l I roo 

MW n h n+M J 2M\n n h n+M 







M!(n - 1)! ~ V n 
where IS* 2 ™" 1 ) = 1ix n /{n — 1)! denotes the volume of the In — 1 dimensional sphere. 
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Turning to I2 in (4.5) we make two changes of variables, z = r9, then r 2 = t, so that 

r 2M+2m+2n-l / 2\ 

|c2n-l| r v^-'J 
t M+m+n-l e -t/h dt 



2Mh n h n+M 



't£[l±£] 

The last integral is very close to the integral appearing in the following lemma which 
will be proved in the appendix: 

Lemma 4.3. Let p > 0, 5 < 1. We define 

J (p, 5) := [ t p e- pt dt. 
Jt(£[l-5,l+8] 

Then 

W)<^exp(-p(l + £)). (4.8) 

To apply this lemma to the last integral in (4.7) we make the change of variable 
t/h = (M + m + n — l)s. To assure that the interval of integration does not change 
much, we claim that under assumptions of Lemma 4.1 we have, 

[1 ± e/2] c — — 1 -[1 ±e] = -[1 ± e]. (4.9) 

h[M + m + n — 1) a 

Indeed, (4.2) implies the following inequalities: 

l-~>^(l-£), 1 + | + ( 4 - 10 ) 

la I a 

The first one is straightforward, since it is equivalent to 2a — 2 > {a — 2)e, and a — 2 < 0. 

The second inequality in (4.10) is equivalent to (2cr — 2)/(2 — a) < e, so that in view 

of (4.2) we need to check that (2a — 2)/(2 — a) < A(a — 1) which follows from the 

assumption a < 3/2. 

Returning to (4.7) we have 

t n+m+M-l e -t/h dt < ^ M + m + n _ i^M+m+n f ^ 
't4[l±e\ Js<£[l±e/2] 

Applying Lemma 4.3 gives 

f £n+m+M-i e -t/h d £ < [h{M + m + n — l)] M+m+n c _( M+n+m _ 1 )( 1+£ 2 /16 - ) 
Jt$[i±e/2] (M + m + n- l)e 2 



[h{M + m + nj\ (M+n+m)(1+£2/16) 

£ 2 



M+m+n 



Hence 



I c2n— II —M—m-n 1 
h < {h{M + m + n)r+m+ »^_Ji_ J_ e -M.V« 

p—M—n—m 1 

< h m (M + m + n)^^——-e-^ 



(4.1i; 
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To simplify the upper bound in (4.11) we first use Stirling's formula to obtain (with a 
small irrelevant loss since k k < k\e k /\fk) 

(M + m + n ) M+m+n < (M + m + n)\e M+m+n . 

Thus the bound in (4.11) can be replaced by 

/ 2 < h m (M + m + n)\ ^ e _ MeVl( . < hm{M + m + n)m+n l_ e -Meyx e _ 

Combining this with the bound (4.6), and applying Lemma 4.2 gives 

\\K\\ < h rn (M + n) n (M + m + n )™+™ J_ e - h ~ 1/3 / 16 

e 2 

< h m (M + n + m) 2n+m \e- h ~ 1/3/w . 

e 2 

This completes the proof of Lemma 4.1. □ 



5. Proof of Theorem 2 

We now combine the basic inequality (3.15) with the estimate on E £ given in Lemma 
4.1. We split (3.15) into four terms: 

(i) Aq = X (/) (1 — 2e) m which is the leading term; 

(ii) A\ = X (/) E e (h, m + N, 0) (1 — 2e) m decreases exponentially to as ft, — >■ 0; 
(Hi) Ai = A (/) — (nm 2 hy (1 + 2e) m ~^ will be estimated by noting that it is 

3=1 J ' 

dominated by its first term; 



[IV 



m j 

A3 = A (/) — (nm 2 h,y E e (h,m + N,m — j) will require more care but de- 



3=0 J 



creases exponentially to as h — ^ 0. 

We want to optimize the parameters h, M, e as functions of the order of /, m, and 
the dimension n. We aim to show that Aq 3> Ai, A2, A3, using Lemma 4.1. For this 
we need to check that the assumption (4.2) is satisfied. 

The basic strategy is outlined as follows 

• (4.2) is satisfied if for all < j < m, h' 1 ~ N + 2m + n — j and h(N + 2m + 
n — j) < 1. Thus we need h~ l ~ N ^> m,n. 

• Aq > A\\ we want to apply Lemma 4.1 and thus we need e 2 /h > —nlog(h); 

• Aq > A2. for this to hold A2 has to be greater than the first term of the sum 
in A 2 , nm 2 (l + 2e) m ~ l h; thus the term (1 + 2e) m has to remain bounded as 
m — > 00: we need e < 1/m. 

• Aq > A 3 : the term A Q has - at least - to be greater than the first term of the 
sum in A 3 ; thus we need to have Aq > E e (h,N + m, m); using Lemma 4.1, 
this holds when e 2 /h > — (n + m) \og{h). 
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We define e as e = h a , where a will be determined. From the considerations above we 
need 

h 2a-l >( n + m ^ log and h a <1 jm. 

To express this as one condition, we demand a = 1 — 2a, that is, a = 1/3. This leads 
to the necessary relations: 

e = h~ 1/3 /16, h<(n + my 3 , N = hT 1 . (5.1) 

Application of the estimates on E £ . To use estimates on E e (h, m + N,m — j) for 

< j < m we need the assumption (4.2) to hold. That means that 

3 

1 < h(N + 2m + n-j - 1) < -, 1 > e > A(h(N + 2m - j + n) - 1). (5.2) 

Since N = 1/h, both inequalities are satisfied for all < j < m if they are satisfied 
for j = 0. Recallingl that e = /i _1 / 3 /16 < 1, this in turn follows from 

4/i(2m + n)<£, h(2m + n) < -. (5.3) 

If h < ^(m + n)~ 3 , then 

85 < 5 1 ' 3 



(m + n) 2 (m + n)' 

which implies (5.3). We conclude that (5.2) holds, hence also (4.2), and hence we can 
apply Lemma 4.1 to E e (h, m + N,m — j), < j < m. 

Final estimate on h. We first start by simplifying A\. Lemma 4.1 shows that 
E e (h,m + N,0) <(n + m + N) 2n e- 2 e' Ns ' !/m < t(3//i) 2n+1 e - h " 1/3/16 . 

Thus 

M < A(/)(l - 2e) m (3/h) 2n+l e- h ~ 1/S/w . (5.4) 
To treat A 2 we note that 

A 2 = A (/) ' (1 + 2£) m " J < A (/) (1 + 2e) m (e Wh - l) . 

But since /i < (n + m) -3 , nm 2 h < 1, and thus exp (nm 2 h) — 1 < nm 2 h, and 

A 2 <A(/)(l + 2e) m nm 2 /i. (5.5) 

We finally treat A3. For that, we need the estimate on i? £ (h,m + N,m — j) proved in 
Lemma 4.1: 

E e (h,m + N,m- j) < h m ~ ] {N + n + 2m - ^^-^-V^ 7 ^ 2 / 16 

< (3A^) 2n 5- 2 (3^) m - J e- (m+Ar)£2/16 
<(3//i) 2ri+1 3 m e"^ 1/3 / 16 . 
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Inserting this in the definition of A3, 

^ = A(/)f^ £ (ft,m + ]V im -j). 
3=0 J ' 

this gives 

A, < A if) (3/h) 2n+1 3-e-^ 173 / 16 

j=0 J ' 

<A(f) (3/h) 2n+1 3 m e"^ 1/3 / 16 . 
Here we used again nm 2 h < 1. Thus we get: 

^3 <3 m (3//i) 2n+1 e-^ 1/3 / 16 . (5.6) 

We recall that we are looking for /i such that for h < h , 

X(f){l-2e) m >A 1 +A 2 + A 3 (5.7) 
is satisfied. In view of (5.4), (5.5), (5.6), to obtain (5.7) it is sufficient to have 

A(/)(l - 2e) m > 3A(/)(1 - 2e) m (3/h) 2n+1 e - h ~ 1/Z/1 \ (5.8) 

A(/)(l - 2e) rn > 3A(/)(1 + 2e) rn nm 2 h, (5.9) 

A(/)(l - 2e) m > 3A(/)3 m (3/h) 2n+1 e - h ~ llz l^. (5.10) 
Since h < 5 = 1/64, e < 1/4 and then (1 - 2e) m > 10 _m ; moreover 

"<a+8sr<('i+-V<i. 



1 — 2e J \ m / 

Thus (5.9), (5.10) can be changed in 

A(/) > 3A(f)nm 2 h, (5.11) 

A(/) > 3A(/)30 m (3/h) 2n+1 e- h ~ ll3/1& . (5.12) 
Since A(/) < A(/), (5.8) and (5.12) are both implied by 

A^j > 3 ■ 30 m (30/h) 2n+1 e"^ 1/3 / 16 . (5.13) 

The logarithmic version of this inequality is 

log yjff) > !og(3) + (m + 2n + 1) log(30) - (2n + 1) \og{h) - /i" 1/3 /16, 

and thus taking h < log (A(/)/A(/)) -3 (m + n)~ 3 log(n)~ 3 assures its validity. Indeed, 

and 

ft < (nlog(n))" 3 =>■ nlog(/i) > -/T 1/3 . (5.15) 
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The estimate (5.11) is straightforward: we need 

h ~^ n " lm "- (5 - 16) 



Let us chose 

h < min 



wr log (w)) ' ( m +^ 3lo ^)~ 3 - 



Then h satisfies the three necessary conditions for Theorem 2 to hold: (5.14), (5.15), 
and (5.16). The bound on iV = 1/h is then given by 

iV > max ( log (jfJ^J (m + n) 3 log 3 n 



which is the same as 



N> j^{m + n) 3 log 3 n. 



Appendix: A non-stationary phase lemma 



We prove Lemma 4.3. Let tp (t) — — log (t) + 1. Then ip is a one to one mapping on 
(0, 1] and on [1, oo). Let us then consider the following integrals: 

1*1—5 poo 

J-(p,5)= e^^dt, J + (p,5)= e^M-^dt. 

JO Jl+5 

The change of variable ip (t) = x gives 



oo / 1 \ —1 



with c~ = ip (1 — 8). Thus we need estimates on p~ 1 (x). But on (0, 1 — 5], we have 
<P (t) < 1 — 5 — log (£). It implies p~ x (x) < e 1 ~ s ~ x . This gives 

/oo e - px 
eX - 1+S _ 1 dx - 

A lower bound for e x ~ 1+s — 1 is 

e ^-l+5 _ 1 > L-l+S _ e -cr\ e x > £ e -l+S+*. 

and hence 

J - (p ,, ) <£^ + ,^._L_i y((1 _ i)e - 1+T 

<-L ((1 -^-r. 



(a.i; 
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The same change of variable applied to J + gives 



00/ 1 \ — 1 



with c + = </? (1 + 5). On (1 + 5, 00), we have <p(t) <t and then (/j" 1 (x) > x. 
Since 5 < 1, 

c+ + < 1_ 

c+ - 1 ~ <p(l + 5) - 1 ~ 5 2 

and thus 

J + (p,S)<^- 2 {(l + 5)e^ s ) p . (A.2) 

Now, 

(l-5)e~ 1+<5 < (l + ^e" 1 " 5 , 5 2 <5, 
and hence the estimates (A.l) and (A.2) give 

J (p, 6) = J_ (p, 5) + J + (p, 5)<^ ((1 + 5) e^Y . 

Also, 

(l + 5)e- 8 <e- s2 /\ 

so that finally 

J(p,5)<-^exp (- p (l + L 
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